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ABSTRACT 

The high redshifts of the most distant known quasars, and the best estimates 
of their black hole masses, require that supermassive black holes (SMBHs) must 
have formed very early in history. Several mechanisms for creating and growing 
these holes have been proposed. Here we present an evaluation of the timescales 
needed for various critical processes in order to discriminate between the proposed 
scenarios. We find in particular that mergers alone are not able to grow the black 
holes at a sufficient rate. Accretion models offer a solution and we use accretion 
timescales to constrain the manner in which the black hole was first formed. This 
analysis implies, but does not require, the action of some unconventional process. 

Subject headings: cosmology: theory — black holes — galaxies: evolution — quasars: 
general 



1. Introduction 

Central supermassive black holes are 
a common feature to galaxies today, but 
which came first, the black hole or the 
galaxy? Conventional thinking would sug- 
gest that the the first generation of stars 
evolved into black holes, which have sub- 
sequently settled to the centers of their 
host galaxies, merged, and accreted gas. 
But this idea, in which central black holes 
form inside pre-existing galaxies, has re- 
cently earned some scrutiny. First, the dis- 
covery of increasingly high redshift quasars 
requires a surprisingly early formation of 
the black holes (see, for example, Fan et 
al. (2001) and Fan et al. (2003)). Second, 



a large quasar sample shows no evidence 
of black holes growing in mass with de- 
creasing redshift (Vestergaard 2002; Diet- 
rich et al. 2002). So we are left to consider 
the possibility that either the central black 
holes formed before their host galaxies, or 
they grew to maturity very quickly within 
them. Either way, they have grown little 
since the quasar epoch. 

The most distant known quasar lies at 
z = 6.41, with a central black hole of 
mass M. = 3 x 10 9 M (Willott, McLure, 
& Jarvis 2003). In the ACDM cosmol- 
ogy observed by WMAP (Bennett et al. 
2003), with n A = 0.73, fi m = 0.27, and 
H = 71 km s -1 Mpc -1 , this redshift cor- 



responds to a time when the universe was 
only 880 Myr old. For the present work, we 
will take this as the time to beat: 3 billion 
solar masses in 880 million years. 

In the past year, two separate HST 
studies have cited evidence for intermedi- 
ate mass black holes (IMBHs) in the cen- 
ters of globular clusters: a 4000 M Q hole 
in M15 (van der Marel et al. 2002), and a 
20 000 M Q hole in Andromeda's Gl cluster 
(Gebhardt, Rich, & Ho 2002). This is the 
lastest and strongest evidence for IMBHs, 
but there is additional evidence, and good 
theoretical motivation as well; see van der 
Marel (2003) for a comprehensive review. 
IMBHs are widely believed to be a neces- 
sary growth stage for SMBHs. In section 2 
of this paper, we will review the major pro- 
posed routes to the formation of a SMBH, 
all of which include an IMBH phase, from 
which time the question is simply one of 
growth. 

We start in Section 2 with a flowchart of 
avenues leading to the creation of a SMBH. 
In Section 3, we examine the timescales for 
each needed process. We conclude in Sec- 
tion 4 by discussing how realistic each av- 
enue is in light of these timescales. 

2. Flowchart 

There are essentially four proposed fam- 
ilies of models leading to the formation of 
IMBHs, and two or three ways to grow 
them. These approaches are depicted in 
figure 1 and discussed in turn below. 

I. The black holes may be primordial, 
in which case they formed from primordial 
density variations before big bang nucle- 
osynthesis. 

Primordial black holes (PBHs) have 



been studied extensively, and the most di- 
rect mechanism for their creation is the 
collapse of gaussian random density fluc- 
tuations (Carr & Hawking 1974). These 
holes come from horizon scale (or smaller) 
modes, and therefore their masses are de- 
termined by their time of formation. In 
the radiation dominated early universe, 

M. ~ 10 5 M . (1) 

But in order to preserve the successful 
BBN prediction of light element abun- 
dances, there must be no significant rate 
of PBH formation once nucleosynthesis be- 
gins, and therefore the PBHs are capped 
at intermediate mass. In addition, Carr, 
Gilbert, & Lidsey (1994) have pointed out 
that, given a small scalar spectral index — 
n s = 0.93 ± 0.03 was recently observed in 
the CMB (Bennett et al. 2003) - PBHs 
from density inhomogeneities should only 
have formed in quantities too small to be 
of interest. 

A more promising, and perhaps in- 
evitable mechanism for forming PBHs 
also exists, in which the collapse is trig- 
gered by "color holes" at the quark-hadron 
phase transition (Crawford & Schramm 
1982). However, because this occurred at 
~ 10~ 6 s, these PBHs would be smaller 
than ~ 1 M Q by eq. 1, and would remain 
as collisionless dark matter today, rather 
than collecting into larger black holes. (In- 
terestingly, Hawkins (1996) shows evidence 
for such PBHs in the microlensing of dis- 
tant quasars, in numbers comparable to 
that needed to account for dark matter.) 

2. Normal population III stars formed 
at ~ 100 M , evolved to black holes, and 
merged at the center of their small dark 
matter halos. 
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Fig. 1. — The various ways to build the earliest supermassive black holes in the universe. 
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This is perhaps the least exotic way 
to create IMBHs, and at this point there 
is very strong support for the process in 
the form of numerical simulations of struc- 
ture formation (Abel, Bryan, & Norman 
2000; Yoshida et al. 2003). These sim- 
ulations include the relevant atomic and 
molecular processes in the first gas clouds, 
particularly cooling by rotation of molecu- 
lar hydrogen, superimposed on CDM halo 
evolution. They find that < 10 stars of 
~ 100 M form in halos of ~ 10 6 M , en- 
gaging < 1% of the system's baryonic mat- 
ter. 

Because the cooling of population III 
stars is hindered by the lack of metals, 
these first stars would be more massive 
than those allowed by fragmented star for- 
mation today. Heger & Woosley (2002) 
have shown that these massive stars will 
evolve into black holes containing a signifi- 
cant fraction of the star's initial mass (this 
fraction varies with the stellar mass, and is 
of order 50%), unless the stars are in the 
mass range 140 M < M < 260 M , in 
which case they are completely disrupted 
when they go supernova. Given the small 
initial size of such a cosmologically young 
halo (~ 10 pc), the holes sink to the center 
and presumably merge into a single IMBH 
there. 

3. Supermassive stars may have been 
the first baryonic objects to form. If 
so, they would have evolved rapidly into 
IMBHs. 

Supermassive stars (SMSs), as a class of 
objects, span from 10 3 to 10 s M , although 
the first generation of them would reside 
at the lower end of that range (Shapiro 
& Teukolsky (1983) provides a comprehen- 
sive introduction on the topic). A 10 3 M 



SMS has a lifetime of < 10 6 years (SMS 
lifetimes range down to 10 years), at the 
end of which, it undergoes relativistic col- 
lapse to a black hole in a matter of seconds 
(New & Shapiro 2001). This collapse sends 
~ 90% of the star's original mass into the 
black hole, with the remaining 10% form- 
ing an accretion disk (Shapiro & Shibata 
2002; Shibata & Shapiro 2002). 

The radius of a 10 3 M SMS, R SM s ~ 
0.1 pc, is around 10% of the virial radius 
for the earliest halos containing gas cool 
enough to collapse. It is therefore reason- 
able to suspect that these stars should ap- 
pear in early halos, especially since pop- 
ulation III stars produce ultraviolet rada- 
tion which splits hydrogen molecules, thus 
destroying the only coolant which could 
collapse smaller stars. The collapsing gas 
might produce enough radiation pressure 
to prevent fragmentation (Baumgarte & 
Shapiro 1999). Another possibility is that 
the first stars formed dense clusters which 
merged into one or more SMSs. However, 
whether or not supermassive stars ended 
the cosmological dark ages remains a mat- 
ter of speculation. 

4. If dense clusters of stars emerged 
from early star forming regions, then the 
cluster stars may have merged with each 
other to create massive stars (say, several 
hundred solar masses), which then evolved 
to massive seed black holes. This idea has 
been mapped out recently by Ebisuzaki et 
al. (2001), and we will not repeat it here. 
However, for the sake of timescales, which 
are the focus of this paper, this scenario 
is nearly the same as the population III 
star scenario discussed above, because the 
merging occurs during the lifetime of the 
cluster's massive stars. 
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Each of the paths in figure 1 sees a 
phase with an IMBH imbedded in a halo, 
at which point accretion and halo mergers 
take over to achieve heavier black holes. 

The merger scenario merges central 
black holes hierarchically along with the 
host halos, growing both together. It is 
clear that merging IMBHs and SMBHs 
are real, interesting phenomena deserving 
study. For example, Milosavljevic & Mer- 
ritt (2001) have simulated the whole pro- 
cess, and shown that coalescing SMBHs 
remove the steep cusps that CDM simu- 
lations generate, thus matching observed 
galaxy rotation curves. But it remained to 
be shown here whether or not mergers are 
prevalent enough to sufficiently increase 
black hole mass. 

Gas accretion is another major way to 
grow IMBHs into SMBHs. The belief that 
quasars are powered by accretion at the 
Eddington luminosity gives this idea cred- 
ibility. But over what fraction of cosmic 
history did these holes accrete? Some 
authors suppose that Eddington limited 
accretion is continuous since the black 
hole first formed (Haiman & Loeb 2001), 
and thus its mass grows rapidly. Others 
have proposed scenarios in which accre- 
tion is negligible over the history of the BH 
(Madau & Rees 2001), such that 150 M 
holes that formed in the first episode of 
star formation still populate our Galactic 
bulge today. 

Another possibility is that the BH feeds 
on dark matter, in which case the halo 
and the IMBH evolve together. MacMillan 
& Henriksen (2002) proposed a model in 
which the halo and hole grow self-similarly, 
with both analytical and numerical results 
(MacMillan & Henriksen 2003) support- 



ing their idea. This model is particu- 
larly interesting because it also simulta- 
neously explains both the spectral index 
for the galaxy power spectrum, n m —2, 
and the M. oc o"b ulgc relationship observed 
to be universal in galaxy centers (Fer- 
rarese & Merritt 2000; Gebhardt et al. 
2000; Tremaine et al. 2002). 

All of these methods certainly con- 
tribute to the growth of massive black 
holes. In the next section we will use 
timescales to help constrain which meth- 
ods could realistically dominate the pro- 
cess of forming SMBHs in less than a gi- 
gayear. 

3. Timescales 

First, we examine the evolution to the 
IMBH phase, as depicted in Figure 1. If 
the BH is primordial, then it can be up 
to 10 5 M Q in the very early universe. Its 
growth by merger or accretion then awaits 
the formation of matter clumps to sur- 
round the PBH. These could begin to 
arrive as early as photon decoupling, so 
timescales alone will not strongly constrain 
this possibility unless the PBH in question 
started at < 1 M Q , as it would if it formed 
during the quark-hadron phase transition. 
In this case the BH seed has a very long 
way to grow, which is barely possible by 
accretion, as will be shown in subsection 
3.3. 

The simulations of population III star 
formation of Yoshida et al. (2003) pro- 
duce some very useful information. The 
earliest simulated star forming region oc- 
curs at z — 32, with stars forming at 
10 -3 M Q yr _1 per comoving Mpc 3 , increas- 
ing to 0.1 M yr _1 Mpc -3 at z = 7, and 
holding constant thereafter (Hernquist & 
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Springel 2003). 

Molecular hydrogen coolant forms in 
these clouds in ~ 30 Myrs. By z = 23, 
the gas pressure has diminished and gas 
clumps begin to contract. Halos are ~ 
5 x 10 5 to 2 x 10 6 M by this time, and the 
largest virialized star forming clouds are 
10 5 M . The gas takes an approximately 
isothermal density profile, with constant 
density inside a core radius of ~ 10 pc. 

The virial temperature is a few thou- 
sand Kelvins. A simple application of 
the virial theorem gives a virial radius of 
45 pc, and the corresponding circular or- 
bital speed is 10 km s -1 . 

This, then, is the starting point for pop- 
ulation III stars leading to IMBHs, with 
star formation beginning in earnest at z = 
23. To convert redshifts to times in a 
ACDM cosmology with f2 = 1 (Peacock 
1999), 

2 Er „ 1 sinh- 1 v /|q n -l|(l + ^)-3^ 

(2) 

and with WMAP values inserted, this is 

t age = 10.8 sinlr 1 (V2.7(l + z)- 3 ) Gyr . 

(3) 

So gas collapses into the first stars at t age ~ 
0.15 Gyr. 

These stars collapse in the free fall time, 

which is about 30 Myr at the densities ob- 
served in simulations by Abel, Bryan, & 
Norman (2000). They evolve to black holes 
(or disruptive supernovae) in a few million 
years, leaving ~ 0.7 Gyr remaining to grow 
into a quasar power source. 



Using the dynamical friction formula 
from Chandrasekhar (1943), quoted here 
from Binney & Tremaine (1987), we can 
find the time needed for these population 
III remnant BHs to sink to the center of 
their host halo. Treating the BHs as point 
particles, 

1.17 vrl 

where In A is called the Coulomb loga- 
rithm, which characterizes the distance 
range perpendicular to the path of mo- 
tion over which the frictional encounters 
occur (for a recent study of this param- 
eter, see Spinnato, Fellhauer, & Portegies 
Zwart (2003)). Using the values mentioned 
above for the parameters in this equation, 
and choosing A m 20 (this choice is com- 
mon and doesn't affect the outcome much 
because of the logarithm), the dynami- 
cal friction timescale relevant here is at 
least ~ 0.1 Gyr, bringing the total down 
to ~ 0.6 Gyr available for BH growth. The 
coalescence time of these holes is not well 
known, and could be long. 

If the star formation yields one or more 
supermassive stars, rather than multiple 
normal population III stars, then the evo- 
lution to a BH is faster than for normal 
stars (Shapiro & Teukolsky 1983). The 
dynamical friction time to merge smaller 
stellar BHs is also avoided with SMSs, so 
in the SMS case, we have ~ 0.7 Gyr re- 
maining to grow a SMBH. 

In the case of an early-forming dense 
star cluster in which runaway collisions 
merge the stars, the timing remains essen- 
tially the same. One waits until the clus- 
ter stars form, merge, and evolve to one or 
more BHs. The time available to become 
supermassive is still at most 0.7 Gyr. 
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3.1. Dynamical Friction 

We consider the case of a satellite sub- 
halo halo falling into a larger halo. Each 
halo is modeled as a singular isothermal 
sphere, so that p oc r~ 2 and orbital veloc- 
ities within a halo are constant. We seek 
a formulation that tracks the friction until 
the satellite's central BH reaches the new 
center. 

Smaller halos formed earlier in the CDM 
hierarchy, when the universe was denser, 
and therefore they are more tightly bound. 
So as a subhalo sinks, we assume it holds 
itself together, until it gets close to the cen- 
ter of the larger halo. There, which halo an 
orbiting particle belongs to becomes am- 
biguous. This occurs roughly when a par- 
ticle's distance from its original host sub- 
halo is on order of the separation between 
the two halo centers. This separation dis- 
tance is r, a radial coordinate extending 
from the center of the larger halo. 

The internal speed within an isothermal 
halo is 



GM 
~R~ 



(6) 



The square root in this equation is help- 
ful because it reduces the dependence that 
velocities have on a halo size R, which is 
generally not well known. 

The mass which can be associated with 
the sinking satellite is then 

2 

M s ~M r + M. = ^ + M., (7) 

where M r is the mass of the satellite halo 
residing inside a distance equal to r around 
it, M. is the satellite subhalo's central BH 
mass, and t> s is the orbital velocity inside 
the satellite. Choosing to use r in this way 
assumes that the isothermal satellite halo 



extends to distances larger than some ini- 
tial infall radius: R s > r^. 

The angular momentum of the satellite 
in the larger halo is L = f x M s v, where 
v is the orbital velocity within the larger 
halo. Dynamical friction induces a torque 
r = fxF = fx M s (dv/dt). 

Following the standard derivation pre- 
sented in Binney & Tremaine (1987), the 
statement of r = dL/dt leads to the equa- 
tion 

dr GM S , . . 

r— = -0.428 In A . 8 

dt v w 

We next define 3 constants as follows: 
Ci = (0.428/?;) In A, C 2 = dGM., and 
C 3 = —C±Vg. Inserting the satellite mass 
from eq. 7, we have 



„ _i „ 
Tt+ C,r ' = £7,. 



(9) 



Eq. 9 is a first-order, linear, ordinary 
differential equation, but the closed-form 
solution is prohibitively long and compli- 
cated due to the r" 1 factor in the second 
term. A numerical solution readily pa- 
rameterizes the sinking orbit r(t), and a 
root can be found for the time tdf at which 
r — > 0. One needs only to specify a start- 
ing radius for the satellite halo, r\ = r(0). 
An example orbit decay curve computed in 
this way is given in figure 2. 

This dynamical friction model approx- 
imately doubles the time spent spiralling 
inward, relative to the Chandrasekhar for- 
mula for an infalling point particle, eq. 
5. Colpi, Mayer, & Governato (1999) 
found the same result using both A^-body 
simulations and the theory of linear re- 
sponse. A subsequent paper by Taffoni et 
al. (2003) gives an analytical approxima- 
tion for NFW halos (Navarro, Frenk, & 
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Fig. 2. — Dynamical friction infall profile 
as calculated in section 3.1. In this figure 
we have exaggerated the curve at the end 
for clarity by using an unrealistically large 
BH mass in the satellite subhalo. 

White 1996), also finding that the decay 
time is increased by a factor of a few if the 
satellite has less than ~ 10% of the larger 
halo's mass. In the mass range of interest 
for the present work, all these models es- 
sentially agree; however, since we wish to 
follow length and mass scales ranging down 
to a central BH, we have chosen to use the 
method described above. 

3.2. Mergers 

We first consider the viability of black 
hole merging as a means to grow SMBHs. 
As halos merge over cosmic history, their 
central BHs sink to the new center by dy- 
namical friction, and subsequently coalesce 
into a larger BH. 

The usual way to track mergers is by 
iV-body simulation. We begin with such a 
simulation, and calculate the total mass of 
subhalos that can merge into one large halo 
between the epoch of the first massive BHs 
and the early quasar epoch. We subject 
this calculation to the constraint that the 
merging subhalos must sink to the center in 



the time allowed between the two epochs. 

This can become a very complicated 
calculation, because the halo population 
is continuously changing over a long pe- 
riod of cosmic time. To make the prob- 
lem tractable, we make a series of major 
approximations. We describe them along 
the way, and discuss them afterward. We 
demonstrate first that with our assump- 
tions, mergers are not sufficient to grow 
SMBHs, regardless of how the BHs origi- 
nally formed. Second, we show that a more 
accurate calculation to replace our simpli- 
fying assumptions would only strengthen 
that conclusion. 

The argument begins with the numer- 
ical simulations of Volonteri, Haardt, & 
Madau (2003). They follow a 4 x 10 13 M 
halo today back in time to all its subha- 
los, presenting halo mass spectra at mul- 
tiple redshifts. For our purposes, the re- 
sults they present at z = 5 and z = 20 are 
most important, approximating the time 
of the first quasars (z = 6.41) and the first 
IMBHs in halos (after z = 23). We'll ad- 
dress the impact of these approximations 
in redshift and total mass after the calcu- 
lation. 

The halo functions are nearly power law 
in the mass range of interest. For M in M Q 
units, at z = 5 



dN 
dM 

and at z = 20 

dN 



10 14 AT* , 



dM 



= 1.6 x 10 16 M 



-2.9 



(10) 



(11) 



Consider the largest halo in existence 
as the period of interest begins. Imagine 
(conservatively) that all other subhalos are 
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available to merge into this one, and (also 
conservatively) that all of these subhalos 
contain central BHs of 0.1% of the par- 
ent galaxy's luminous mass, as observed 
in large galaxies today (Ferrarese & Mer- 
ritt 2002). According to WMAP, the ratio 
of baryonic matter to total matter is 0.17 
(Bennett et al. 2003), so conservatively we 
may say that a central BH has 0.02% of the 
host subhalo's mass. Then every subhalo 
massive enough to dynamically sink to the 
center in the available time does sink, lead- 
ing to a more massive combined halo. 

If all the IMBHs from all the combined 
subhalos coalesce quickly (again, highly 
conservative), then the new SMBH mass is 
0.02% of the total mass of all subhalos with 
mass M > M s min , the minimum satellite 
mass needed for dynamical friction to sink 
a subhalo by the early quasar era. 

However, during the merging process, 
the largest halo is gaining mass and the 
number of massive halos is increasing. To 
evade this problem, we suppose that the 
halo mass spectrum at the end of the merg- 
ing process (at z — 5) applies through- 
out the merging process. Our calculations 
show that in most cases, this assumtion 
conservatively exaggerates the final SMBH 
mass by more than an order of magnitude. 

Using eq. 10, the total mass of the 
quasar's SMBH is 

POD 

M. = 0.0002 / 10 14 M"5 dM . (12) 

The lower limit, M Sjmin , is determined by 
dynamical friction as described in section 
3.1. We follow the infalling halo mass, 
and transition to the (smaller) infalling BH 
mass as the two holes close in, achieving re- 
sults closely consistent with other models 



that include mass loss but not a central BH 
(Colpi, Mayer, & Governato 1999). 

In the dynamical friction calculation, 
three parameters must be chosen. We take 
the internal velocity of the larger halo to 
be v 100 km s _1 , noting that for a fixed 
dynamical sinking time (fixed by the time 
available for BH growth), the final SMBH 
mass obeys M, oc t> -1 / 3 . So choosing a 
larger velocity would only very slightly re- 
duce the final SMBH mass by reducing the 
number of subhalos which sink in the al- 
lowed time. 

The other needed quantities are radii: 
Ti, the initial radius from which the satel- 
lite halo begins its descent, and R s , the 
size of the satellite subhalo. The former is 
roughly the impact parameter of the col- 
liding subhalos, and the latter indicates 
the subhalo's concentration factor. Both 
radii will clearly vary from one encounter 
to the next, so we choose several repre- 
sentative values. The results demonstrate 
that any reasonable (in fact, conservative) 
choice of these numbers leads to the same 
conclusion, that merging black holes take 
too long. 

As shown earlier, the timescale between 
the first IMBH and a quasar's SMBH de- 
pends on how the BH formed. A primor- 
dial BH has the age of the universe at the 
time of the most distant quasar's emission 
known today, or 0.88 Gyr. Population III 
stars and supermassive stars have 0.7 Gyr. 
If the population III stars merged within 
their original halo first, then they have 
0.6 Gyr to grow. 

We fix the sinking time at the appropri- 
ate value and compute the total BH mass 
that can be assembled by halo mergers, ne- 
glecting the time needed for coalescence. 
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Table 1: SMBH masses achieved through 
halo mergers by the early quasar epoch. 
Rs is a characteristic size for the satellite 
subhalos, and rj is the radius from which 
they are assumed to begin their infall. t is 
the time allowed for growth. 

Our results are listed in Table 1, which 
lists a maximum SMBH mass for various 
choices of i? s and r^. No choice brings M. 
up near 3 x 10 9 M , the value needed to 
explain the most distant known quasars. 

One is led to conclude that mergers 
alone are not sufficient to grow the largest 
SMBHs in the oldest quasars. Note that al- 
though the argument above relies on many 
simplifying assumptions, all but one tend 
to overestimate the resulting SMBH mass, 
so the SMBH masses in table 1 should be 
considered extreme upper limits. Our as- 
sumptions are summarized as follows. 

1. All halos contain central BHs at 
0.02% of the halo mass (this may seriously 



overestimate the number of BHs). 

2. All relevant merging is completed by 
z — 5, instead of today, so the number of 
infalling BHs is again exaggerated. 

3. In the ACDM world, going from red- 
shift 20 to 5 allows over 1 Gyr for merg- 
ing, so the number of more massive subha- 
los available to merge is artificially large. 
This overestimates the final SMBH mass, 
because larger infalling subhalos will dy- 
namically sink faster. 

4. Infall radii begin at < R s , which 
for some mergers represents a head start of 
multiple orders of magnitude. 

5. The choice of an isothermal halo 
model probably exaggerates the central 
density, artificially speeding up the action 
of dynamical friction near the center. 

6. SMBH coalescence is assumed to oc- 
cur arbitrarily quickly, but in reality takes 
at least ~ 10 8 years, and probably much 
longer (Milosavljevic & Merritt 2001). 

7. The choice to trace a 4x 10 13 M halo 
today back in time tends to underestimate 
the number of progenitor subhalos. Ex- 
trapolating this mass linearly from Volon- 
teri, Haardt, & Madau (2003), we find that 
increasing it by 100 times would increase 
dN/dM at z = 5 by a factor of < 10 while 
approximately maintaining its slope. 

In this case, the largest halo at z = 5 
would be ~ 10 13 M , so halo speeds should 
be comparable to those in M87 today, the 
giant elliptical at the center of the Virgo 
galaxy cluster. Radial velocities for glob- 
ular clusters around M87 (Cohen 2000) 
range from 1000 to 2000 km s -1 . These 
chane; 1000 km s _1 , allow 

a SMBH to become almost 10 9 M in 
0.7 Gyr from an initial radius r\ = R s = 
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1 kpc, which is a factor of ~ 4 greater than 
our original result. 

So if we attempt to correct for this sev- 
enth assumption, we approach (but do not 
reach) the 3 x 10 9 M hole that quasar ob- 
servations require. But this is only true if 
we take literally the six highly conservative 
assumptions above. 

3.3. Accretion 

Consider first gas accretion, followed by 
feeding on dark matter. 

Gas accretion is probably the simplest 
and most successful way to explain the 
growth of an IMBH into a SMBH, or even 
a large stellar mass BH into a SMBH. 
Haiman & Loeb (2001) show that if the 
black hole accretes at the Eddington limit 
continuously, from birth through detection 
as a quasar, it achieves the needed growth. 
They give the Eddingtion luminosity, 

Le = A.itGM % cii e m v G^} , (13) 

and introduce a radiative efficiency e = 
L/M.c 2 and a fraction of Eddington out- 
put rj = L/Le- The resulting e-folding 
time is 

W = ^ = 4xl0 7 (^)^ 1 yr, (14) 
which equals a 10-folding timescale 

Wo = 9.2xl0 7 (^-)77- 1 yr- (15) 

For a fiducial value of e ps 0.1 and con- 
tinuous Eddington-limited accretion (rj ps 
1), one can calculate the initial BH mass 
needed to grow to 3 x 10 9 M in the alloted 
time. 

If the seed BH is primordial, it might 
accrete steadily from approximately the 



beginning of time. In this case, it has 
0.88 Gyr to grow; that's 9.6 orders of mag- 
nitude in mass by eq. 15. So the ini- 
tial seed mass must have been M seed > 
0.8 M , which, interestingly, could have 
been formed during the cosmological quark- 
hadron phase transition (Jedamzik & 
Niemeyer 1999). 

If the seed BH is a normal population III 
stellar remnant, or a supermassive star's 
remnant, then it had 0.7 Gyr. It's initial 
mass was then M see d > 70 M , which is 
completely plausible for the first genera- 
tion of stars. One star of > 260 M would 
suffice (Heger & Woosley 2002). 

In 0.6 Gyr, the seed BH mass needs 
to be M seed > 900 M . This applies for 
the case where population III stars must 
sink to the center of their parent halos 
and merge before beginning to grow sig- 
nificantly by gas accretion. 

In the case where the seed BH formed 
from the evolution of normal population 
III stars, the mass of the seed star is sug- 
gestively close to what one would expect. 
Thus this route is possible, but only if 
Eddington limited accretion is maintained 
during the entire IMBH growth process. If 
the seed hole evolved from a more massive 
SMS, then there is time to spare, and the 
accretion rate could have dropped for some 
of the growth period. 

The success of this quick exponential 
growth calculation seems to indicate that 
baryon accretion is the most realistic mech- 
anism for IMBH growth. When the BH is 
small, the Eddington mass accretion rate is 
not prohibitively large, and by z = 6.41 we 
have very strong evidence that the quasar 
is accreting at its Eddington rate (Willott, 
McLure, & Jarvis 2003). So it is reasonable 
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to speculate that it has been accreting at 
this rate all along. 

However it is worth noting that the cen- 
tral SMBHs in modern galaxies appear to 
grow by a factor of 10 on timescales of 
> 8 Gyr (Merrifield, Forbes, & Terlevich 
2000). This is very slow compared to 
Eddington limited growth, requiring rj < 
0.01. This dramatic drop in accretion rate 
needs explaining, especially if we are to 
claim that no similar drop has ever hap- 
pened before z ~ 6 in all of cosmic history, 
as the notion of SMBH growth by g 
cretion onto population III stellar corpses 
would mandate. This problem is only alle- 
viated if the first stars were supermassive, 
or the seed BHs were primordial and sig- 
nificantly larger than 1 M . 

Accreting dark matter is another possi- 
bility. MacMillan & Henriksen (2002) have 
proposed an interesting way to account for 
IMBH growth in which the hole grows pro- 
portionally to the dark matter halo as mat- 
ter falls in. As mentioned earlier in section 
2, this assumption leads to the observed 
galaxy power spectrum and the observed 
M. — a relation. 

In numerical simulations (MacMillan & 
Henriksen 2003), dark matter was found to 
form a self-similar region surrounding the 
BH, which does grow in proportion to the 
BH growth at the center. The formation 
of this region proceeds by gravitational in- 
teractions between clumps of dark mat- 
ter ("particles" in the simulations) which 
must have formed earlier in the CDM uni- 
verse. The BH grows whenever the parti- 
cles cross inside the Schwarzschild radius. 
(The BH growth rate is thus understated 
in the sense that i? sc h was used instead of 
the larger capture radius.) 



In the absence of angular momentum, 
the black hole grows rapidly: M. oc t 4 . But 
for realistic departures from pure spheri- 
cally symmetric infall, in which the dark 
matter has angular momentum to help it 
resist falling into the hole, the growth of 
the BH is only linear in time: M. oc t (al- 
though the M. — a relation is recovered 
either way). 

The simulations were not performed in 
standard units, and therefore cannot be 
immediately applied to relevant halo sizes 
and cosmic timescales. Here we simply 
note that the simulated BH growth is slow 
(e.g., linear in time, as compared with ex- 
ponential in time for gas accretion). One 
therefore needs an unusually large IMBH 
to start from, which requires an unusually 
large SMS as its origin. 

4. Conclusions 

In this paper we have surveyed the 
state of knowledge of SMBH formation and 
growth in light of a difficult new data point 
to satisfy: a quasar at z — 6.41 whose 
SMBH was 3 x 10 9 M when the universe 
was only 880 Myr old. We have extended 
the calculations of others on SMBH growth 
via mergers and accretion. 

In the case of SMBH growth by merg- 
ers, we follow the output from a CDM sim- 
ulation in order to track progenitor subha- 
los. We stipulate that a subhalo must sink 
to the center by dynamical friction after a 
merger within the time available for SMBH 
growth. This limits the mass of subhalos 
and therefore BHs which could have coa- 
lesced in time, thus constraining the final 
SMBH mass. We find that in any realis- 
tic case, mergers are incapable of growing 
a BH with the needed speed. 
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Accreting baryonic matter represents a 
viable growth process. For this mecha- 
nism, we impose time limits sensitive to the 
manner in which the seed BH first formed. 
The following possibilities emerge: 

1. A primordial black hole that formed 
during the quark-hadron phase transition 
has accreted at the Eddington rate contin- 
uously ever since baryons have been able 
to cluster until at least the quasar epoch. 

2. A large population III seed star with 
M ps 260 M evolved to a BH, and ac- 
creted at the maximum (Eddington) rate 
continuously between its formation and the 
quasar epoch. 

3. A collection of merged population III 
stellar black holes formed an IMBH with 
M ~ 900 M , and accreted at the Edding- 
ton rate continuously through the quasar 
epoch. 

Each of these options involves seed stars 
at the upper end of their predicted mass 
range. So although the Eddington limited 
gas accretion scenario is adequate to grow 
the holes in the time available, it is not en- 
tirely satisfying, because it requires heavy 
gas infall which is never significantly di- 
minished before the quasar era, despite the 
relative gas suffocation which has evidently 
happened since. One can avoid imposing 
this requirement by starting with a larger 
BH: 

4. A single SMS, or the merged remains 
of a dense cluster of stars, or a PBH formed 
well after the quark-hadron phase transi- 
tion, yielded an IMBH much larger than 
~ 100 M . It accreted through the quasar 
epoch, but did not need to maintain the 
maximum accretion rate during that time. 

One other option, wherein the BH grows 



in proportion to its host halo by eating 
dark matter, is expected to yield slower 
BH growth. As such, normal population 
III stellar seeds would be inadequate, al- 
though a particularly massive SMS seed re- 
mains a possibility. 

We conclude by suggesting that some 
unconventional mechanism is needed to re- 
alistically beat the clock for the early for- 
mation of SMBHs. BH growth through 
mergers is too slow, and growth by accret- 
ing nonstop at the Eddington limit is prob- 
ably too contrived. The remedy appears 
to require one of the following: (a) PBHs 
from unusually large overdensity pockets, 
(b) population III starbursts well in excess 
of those predicted by current simulations, 
or (c) supermassive stars. 

This work was supported by NASA un- 
der the Colorado Space Grant Consortium. 
We are grateful to R. N. Henriksen and J. 
D. MacMillan for helpful correspondence. 
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